Entangled Wavefunctions from Classical Oscillator Amplitudes 
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In the first days of quantum mechanics Dirac pointed out an analogy between the time-dependent 
coefficients of an expansion of the Schrodinger equation and the classical position and momentum 
variables solving Hamilton's equations. Here it is shown that the analogy can be made an equivalence 
in that, in principle, systems of classical oscillators can be constructed whose position and momenta 
variables form time-dependent amplitudes which are identical to the complex quantum amplitudes 
of the coupled wavefunction of an iV-level quantum system with real coupling matrix elements. 
Hence classical motion can reproduce quantum coherence. 
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I. INTRODUCTION 

In their first formulation of quantum mechanics, both 
Schrodinger and Dirac were strongly influenced by con- 
nections to the Hamiltonian formulation of classical me- 
chanics. Indeed, in one of the very first applications of 
Schrodinger's time-dependent equation (TDSE), Dirac 
[l[ indicated a close parallel between the coupled first- 
order set of equations arising from the TDSE and the 
coupled first-order Hamilton equations of classical me- 
chanics. Dirac introduced the time-dependent basis set 
expansion, 



*(*)) =X)c„(*)k* 



(i) 



where the c n are complex co-efficients and | n n ) denotes 
an arbitrary basis. In the TDSE (where, for the moment, 
we put H = 1) this expansion leads to the set of first-order 
coupled equations, 



(2) 



where H nm are the matrix elements of the quantum 
Hamiltonian. Dirac then remarked that by considering 
the co-efflcients q n = c n and p n = ic* n to be canonical 
variables and assuming a 'Hamiltonian function' 



(3) 



the quantum equations are equivalent to the classical 
Hamilton equations 



dU 

9p„ 



Pi, 



dU 



(4) 



Note, however, that this is still fully quantum mechani- 
cal since the matrix elements and the amplitudes appear- 
ing in the 'classical' Hamiltonian arc all complex objects. 



'Electronic address: 
t Electronic address: 



briggs@physik.uni-freiburg.de 



Hence, there is no obvious classical counterpart. Dirac 
also made a transformation to real variables by using am- 
plitude and phase of the quantum co-efficients, i.e. 



(5) 



wire 



re p n 



\c n . Although with this transformation the 
variables are real, the 'classical' Hamiltonian still con- 
tains complex quantities and hence does not obviously 
correspond to any real physical system. 

Forty years later, Strocchi Q approached the question 
of quantum/classical equivalence slightly differently. He 
first formulated classical mechanics in terms of complex 
variables z n = (q n + ip n )/y/2, with p n and q n real, to 
give the Hamiltonian equations Q in the form, 



&H_ 

dz* ' 



(6) 



and its complex conjugate. Here H. is a real classical 
Hamiltonian. Then it is remarked that the TDSE cou- 
pled equations @ arc equivalent to these classical equa- 
tions i.e. c„ = z n if the 'classical' Hamiltonian function 
is taken as the expectation value of the quantum Hamil- 
tonian i.e. 



w = <*(t)|fT|*(f); 



i H nrn c m 

(*)• (7) 



as in the treatment of Dirac. Although ostensibly time- 
dependent and complex, it is easy to prove that % be- 
comes constant and real for Hermitian Hamiltonians. 
Identifying c„ with z n and reverting to the real variables 
(iniPn) one has however , still the apparently complex 
form quoted by Strocchi 



H 
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Of course, for H Hermitian one does achieve a real 'clas- 
sical' Hamiltonian 

T-L = r ^[H mn (q n q m + p n Pm) + 2^s(H nm )q n p m ] (9) 
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and in the special case that all coupling matrix elements 
are real the Hamiltonian 



H 



(q n q m + PnPm), 



(10) 



which is that of coupled real harmonic oscillators. Note 
that the coupling is of a very special form in which there 
is both linear position and momentum off-diagonal cou- 
pling with exactly the same coupling strengths. 

The mapping of the TDSE onto classical-like equa- 
tions up until now has been regarded as something of 
a curiosity, although it has been used as a starting point 
to treat some problems of molecular electronic dynam- 
ics semi-classically 0, H| . Here we show that systems of 
coupled classical oscillators corresponding to the Hamil- 
tonian Eq. (fTOj) can be realised whose position and mo- 
mentum variables reproduce exactly the time-dependent 
coefficients of an expansion of the quantum wavefunc- 
tion. We call these the p&q-coupled oscillators. Each 
oscillator plays the role of a state in Hilbert space and 
the coupling between the oscillators the role of quan- 
tum coupling matrix elements. Unfortunately, for a gen- 
eral quantum system the corresponding exactly equiv- 
alent classical system involves rather complicated cou- 
pling schemes between the oscillators. However, we show 
that in a weak-coupling approximation, simpler sets of 
coupled oscillators in which only the q-coupling in the 
Hamiltonian Eq. (|10[) occurs provide an exceedingly good 
approximation to the exact result. 

In previous papers [H, @, where we studied energy 
transfer along dipole-dipole interacting molecules, we 
called this approximation the 'realistic-coupling approx- 
imation' (RCA) since for this case, either classical or 
quantum, to be realistic it is essential that the oscilla- 
tors and quantum entities largely retain their character 
when coupled i.e. the coupling is weak compared to their 
internal forces. Although in this work we consider a gen- 
eral Hamiltonian where the diagonal elements in a chosen 
basis are not necessarily associated with concrete phys- 
ical entities and thus there is no a priori reason for the 
off-diagonal elements to be small, we will still keep the 
term RCA for this weak coupling approximation. 

It is interesting that, with reference to specific sys- 
tems, several authors have already recognised the similar- 
ity between q-coupled classical oscillators and few-level 
quantum systems. Indeed, already in one of the the first 
quantum treatments of resonant electronic energy trans- 
fer between identical atoms, Frenkel in 1930 Q remarked 
on the essential equivalence of the quantum treatment to 
the earlier 1925 classical treatment of Holtsmark §. In 
particular he showed that the normal mode frequencies 
of the classical coupled equations reduce to the quantum 
eigenvalues in the approximation that we call the RCA. 

In 1977 McKibben § constructed a system of three 
coupled mechanical pendula and showed that, in a weak- 
coupling approximation similar to our RCA, their equa- 
tions of motion become equivalent to the quantum equa- 
tions describing the operation of a spin filter on three 



levels of the hydrogen atom. Later Hcmmer and Pren- 
tiss [To[ used the same 3-pendulum classical system to 
interpret the quite different 3-level atomic problem of 
the stimulated resonance Raman effect and pointed out 
"the strong mathematical similarities between the pen- 
dulum amplitude equations and the resonance Raman 
equations" . More recently, Marx and Glaser [ll[ have de- 
rived an exact correspondence (in the RCA) between the 
dynamics of three isotropically coupled spins and those 
of three coupled classical oscillators. Here the Liouville 
equation for the density matrix and spin correlation func- 
tions derived from it were considered. 

In a rather different application, Jolk et al. [l2| and 
Kovaleva et al. [l3[ have pointed out the similarity of 
the non-crossing behaviour of the eigenfrequencies of a 
pair of dissimilar classical oscillators under varying cou- 
pling to the non-crossing behaviour of the eigenenergies 
of two quantum levels subject to a varying perturbation, 
as in the celebrated Landau-Zener problem. Similarly 
Spreeuw et al. 14| has discussed a classical optical sys- 



tem as analogue to a driven two-level atom. 

In previous communications 0, @ we have shown that 
the time-dependent Schrodinger equation for an aggre- 
gate of N coupled monomers having a single electronic 
transition and the classical equations for coupled elec- 
tric dipoles are equivalent in the RCA. There we demon- 
strated that the quantum coherence in the transfer of 
electronic excitation along a linear chain of monomers 
is reproduced by the transfer of electric dipolc strength 
along the equivalent classical array. The equivalence 
proof was restricted to this special problem of coupled 
classical electric dipoles and its relevance for energy 
transfer in photosynthesis. 

These many and varied applications point to a gen- 
eral property of equivalence of classical oscillators and 
quantum systems in finite Hilbert spaces. However, none 
of the above authors discuss their findings in the con- 
text of Dirac and Strocchi's analysis. Also they all use 
classical oscillators that are q-coupled, i.e. they work 
only in the RCA and do not recognise that it is possible 
to construct classical oscillators that exactly mimic the 
quantum system. In the present paper we show how to 
construct such p&q-coupled oscillators and then clarify 
under what conditions the RCA allows simpler q-coupled 
sets of oscillators to be used as analogues of coupled-state 
quantum systems. 

The structure of the paper is as follows. In Sec. 2.1 we 
derive the classical p&q-coupled equations whose vari- 
ables give directly the coefficients of the time-dependent 
Schrodinger wavefunction. We also examine the deriva- 
tion of the RCA and the meaning of this approximation 
in the context of the Dirac and Strocchi mapping of the 
Schrodinger equation to Hamilton's equations. As spe- 
cific example, in Section 3 we consider excitation transfer 
on a quantum aggregate of coupled monomers and its 
classical equivalent of N coupled pendula. Special atten- 
tion is given to the N = 2 dimer case whose wavefunction 
is considered a fundamental example of quantum entan- 



3 



glement. In this case the quantum (and therefore also the 
p&q-coupled classical) solution can be obtained in simple 
closed from. Then the RCA is tested by comparing the q- 
coupled classical solution with the p&q-coupled solution 
i.e. with the exact quantum solution. Oscillatory energy 
transfer on the dimer is shown to correspond classically 
to the well-known " beating" phenomenon of coupled pen- 
dula from whose motion the entangled wavefunction can 
be extracted. 

The classical analysis is given, as example, for the spe- 
cific case of coupled undamped vertical pendula. How- 
ever it is clear that the analysis also applies to other 
sets of oscillators, particularly classical electrical LC cir- 
cuits, as shown in appendix [5] It is also demonstrated 
in appendix [Cl that the inclusion of a velocity-dependent 
damping term in the classical equations corresponds ex- 
actly to the inclusion of damping via a complex eigen- 
value in the quantum equations, as used for example in 
Ref. 0. 



II. QUANTUM iV-LEVEL PROBLEM AND 
COUPLED CLASSICAL OSCILLATORS 

A. Exact mapping of quantum to classical motion 

First we will show that the first-order Schrodinger 
equations ([2]) map simply onto the classical Hamilton 
equations for a set of coupled classical oscillators. Since 
the classical problem of coupled oscillators is almost al- 
ways formulated using the Newton equations rather than 
the Hamilton equations we also present the correspond- 
ing second-order Newton equations. We start with the 
Strocchi classical Hamiltonian, Eq. (fTU|) . The Hamilton 
equations ^ give 



i Prn Pr\ 



nmQm 



(ii) 



Symbolically, writing q and p as vectors and H as a ma- 
trix the above equations are 



q = Hp 



Hq 



and formally 



Hp 



H 2 c 



(12) 



(13) 



which arc a set of coupled oscillator equations and can 
be solved for q(t) and q(i). The momenta at time t can 
then be calculated from 



p = H 1 q 



(14) 



From the Hamilton Eqs. (|12l) we have 

z = H 2 z (16) 
Similarly the Schrodinger equation ((2|) is written 

it = He (17) 



H 2 c 



(18) 



which is exactly the classical equation (fTSj). Hence the 
p&q-coupled classical equations and the coupled quan- 
tum Schrodinger equations arc identical and in particular 
they have the same eigenvalues. 

To see how to construct a set of real classical oscilla- 
tors, now we write the formal solution in terms of the 
individual oscillator amplitudes. First we split off the 
diagonal term in the Hamiltonian. This quantum en- 
ergy H nn (divided by H) will be denoted by u> n . The 
off-diagonal elements (divided by H) will be denoted by 
V mn ■ Then the Hamilton equations give 



q n = LO n p n + 2J V nm p m 
Pn — ^nQn ^ ^ VnmQm 

Forming z n = \{q n + ip n ) gives 



(19) 



nm^m 



(20) 



which arc identical to the quantum Eqs. ([2]) by construc- 
tion. Taking the time derivative of the q n equation and 
substituting for the p n from the second equation leads to 
the coupled second-order equations 



q n + ui n q n = - 



where we defined 



E 



(21) 



(22) 



This set of classical equations are the p&q-coupled equa- 
tions. Although the coupling elements appear quite com- 
plicated, for a given set of quantum matrix elements V nm 
they can, in principle, be solved to obtain q n (t) and q n (t). 
Then the time-dependent momenta must be calculated 
by matrix inversion as in Eq. (|14[) and the complex z n (t), 
equal to the quantum coefficients c n {t), calculated. 



The set of complex amplitudes, the vector z, is con- 
structed as 



Hamiltonian in RCA 



z=i=(q + l p). 



(15) 



To construct a classical analogue to a given quantum 
system the oscillators must be coupled as in Eq. (|21l) . 
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However, when the Hamiltonian involves only positional 
q couplings, the Hamilton equations simplify to 

q n = LO n Pn 

Pn ^nPn ^ ^ H nm q rn ( ) 

or equivalently 

q n + u? n q ri = - u n V nrn q mi (24) 

ra^n 

which are the q-coupled classical equations. If we 
define real coupling elements K nm = —u)nV nm these 
arc the standard coupled equations, for example of 
a set of linearly-coupled mechanical or capacitatively 
-coupled electrical oscillators. In Appendix A it is 
shown that in the RCA, which corresponds to having 
V nm "C (jJ n ,ui m , Vn, m, the solutions of the simple 
q-coupled equations are a good approximation to the 
more complicated p&q-coupled equations (|21[) 



B. Eigenfunctions and Eigenmodes 

1. Quantum result 

In the basis | ir n ) defined in eq. (UJ the Hamiltonian 
defined in is not diagonal in the general case. To find 
the eigenfunctions and eigencnergies we solve 

H\i> k ) = E k \iP k ) (25) 

The eigenstates | ip k ) can be expressed in the original 
basis via 

N 

\^k) = J2 B kn\^n), (26) 
n=l 

Denoting the diagonal elements H nn = ( ir n \H\ ir n ) = e n 
and the off diagonal elements by H nm = ( ir n \H | 7r m ) = 
Vnm the co-cfficicnts B kn can be obtained from the cou- 
pled set of equations 

(E k — e n )Bkn = ^ VnmBkm (27) 
m 

A general time-dependent wavefunction (coherent 
wavepacket) can be expanded either in the basis 1 7r„ ) 
or in the eigenbasis \ipk), i-e. 

i*w) = E c «wi^) = E^wi^) ( 28 ) 

n k 

Since the eigenbasis diagonalises the Hamiltonian, the 
coefficients b k (t) are given simply by, 

b k (t)=A k exp(-lE k t] , (29) 



where the time-independent complex coefficients A k are 
decided by the initial conditions. Then using Eq. (|26p in 
Eq. ([28]) one has the amplitude of state n, 

Cn{t) = J2 AkBkn ex P \\ Ekt ) ■ (30) 

k ^ ' 



2. Exact mapping 

The Newton equations of motion for the exactly equiv- 
alent set of classical oscillators are the coupled equations 
Eqs. (|2"Tj) . The eigenfrequencies of these classical equa- 
tions are identical to the E k of the diagonalised quantum 
problem. 

Since the eigenvalues and therefore the eigenfunctions 
of the quantum and classical systems are identical, we see 
that the classical time dependence can be used to con- 
struct the quantum wavefunction. The normal modes 
arising from the diagonalization of Eqs. (f2~Tj) are of the 
form q k = (3 k cos(Q, k t + a k ) where /3 k and a k are real con- 
stants. From these normal modes we derive the velocities 

q k (t) = -n k /3 k sm(n k t + a k ), (31) 

Then the complex classical amplitudes are z k (t) = 
-^{qk(t) + ip k (t)) = -^(qk{t) + (i/n k )q k {t)) to obtain 
the general solution as a sum of normal modes, 

*.(*) = E BknAk expH(Ofct)), (32) 
k 

where wc have to set the t — initial conditions so that 
A k = {(3 k /^/2) exp(zafc) and put tt k = E k /h. Comparing 
this result with Eq. ([30]) we see that the classical am- 
plitudes z n (t) reproduce the quantum co-efficients c n (t) 
and the classical motion can be used to reconstruct the 
quantum entangled wavefunction of Eq. (1281) . 



3. Coupled equations in RCA 

As we show in Appendix A, in the RCA the Eqs. (|2"Tj) 
reduce to become identical in form to the simpler linearly 
g-coupled classical equations (f2"4")l , which we write, 

q n + uj 2 n x n = - ^ K nm q m . (33) 

Now we show that, again in RCA, these equations have 
eigenfrequencies which closely approximate those of the 
quantum problem and therefore also the eigenfrequen- 
cies of Eqs. (|2~Tj) . Substitution of the special eigenmode 
solution q n = C kn cos (tt k t) gives 

(fig - u 2 n )C kn = £ K nm C km . (34) 

m 
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This equation can be written 

(fife - U n )C kn = (fife + W„) _1 ^ K nmC km - (35) 

rn 

Now we make a second time the realistic coupling approx- 
imation (RCA), by considering that the spread of eigen- 
frequencies (bandwidth of the dispersion relation) result- 
ing from diagonalisation of Eq. (f34|) is small compared to 
the mean natural frequency to = uj n of the oscillators. 
Similarly, for non-identical oscillators, the width of the 
lu„ distribution must be small compared to lu. Then on 
the r.h.s. of Eq. (|34[) wc can approximate fife and tu n by 
lu to give, 

(tt u -u n )C kn = Y,^C km . (36) 

m 

Further, if we make the identification E k = h£l k , e„ = 
huj n and V nm /h = K nm /(2u), Eqs. (|2"T1) and are 
identical and we can put in RCA 

C kn = B kn V k, n. (37) 

Note that the RCA is synonymous with the condition 
that the couplings K nm are small compared to lu and 
cquivalcntly in the quantum case that all V nm , are small 
compared to the mean energy e of the various eigenvalues 
e n of H . Note also that this in no way implies that 
perturbation theory must be applicable since, as in the 
examples given below, the analysis applies to the case 
where all e„ are degenerate, where perturbation theory 
is invalid. 

III. TWO SIMPLE EXAMPLES 

To illustrate the construction of classical oscillator 
systems which can reproduce entangled wavefunctions, 
both exactly and approximately in RCA, we consider 
the quantum problem of N identical two-level systems. 
Specifically we consider the dynamics when there is ex- 
actly one excitation present and take 1 7r„ ) as the state 
in which monomer n is excited and all other monomers 
are in the ground state. This is the exciton model stud- 
ied originally by Frenkel. Such model quantum systems 
are of fundamental importance in several areas, for ex- 
ample in quantum computing and in the modelling of 
the photosynthetic unit, in addition to applications men- 
tioned in the introduction. Indeed, the quantum dimer of 
N = 2 coupled two- level systems (two qubits), is viewed 
as having a wavefunction which is the simplest example 
of quantum entanglement. This is the case we discuss 
first. 

A. The quantum dimer problem and two classical 
coupled oscillators 

To make the problem concrete we will think of a dimer 
composed of two identical atoms or molecules each hav- 



ing only a ground and one excited state. Note, that we 
restrict to the subspace where exactly one excitation is 
present. Classically this corresponds to two coupled os- 
cillators. 



1. Quantum result 

The dimer has + and — cigenstates of the form, 

|V±> = ^=(K>±k 2 » (38) 

with eigenenergies e± — e ± V, where e is the monomer 
transition energy and V = V12 = V21 is real. The 
state 1 7Ti ) has monomer 1 excited and monomer 2 in 
the ground state and correspondingly for 1 772 ) - 
We expand a solution of the TDSE as 

|tt(t)>=o+(f)|^ + >+a_(f)|V>-) ( 39 ) 

The two eigenstates propagate independently in time ac- 
cording to, 

a±(t) =A± exp[-(i/h)e±t] = A± exp[-(i/h)(e± V)t] 

(40) 

Excitation transfer is described by the initial condition 
I ^(0) ) = 1 7Ti ) which leads to the time-dependence, 

|*(t))= Cl (t)|7ri> + C2(t)|7r 2 ) (41) 

with, C1/2 = h [ CX P ( — ie+t/K) ± exp (— ie-t/H)] . These 
expressions can be simplified further to give 

ci(t) =exp[-(i/K)et]cos[Vt/h] (42) 
c 2 (t) = -iexp[-(i/h)eb] sm[Vt/h] 

which are the exact quantum solutions and describe a pe- 
riodic transfer of excitation between the two monomers. 
Note that in forming the density matrix of coefficients 
c*cj the pure phase factor exp[—(i/h)et] disappears. 

2. Exact mapping 

The mapped Hamilton equations, with uj = e/H and 
V = V/H, for the case N = 2 are 

qi = topi + Vp 2 pi = -ojqi - Vq 2 

(43) 

q 2 = Lop 2 + Vpi p 2 = -tuq 2 - Vqi 

to give the coupled oscillator equations 

q\ /2 + (u 2 + V 2 )q 1/2 = -2u;Vq 2/1 (44) 

In the usual way these symmetric equations can be di- 
agonalised by the transformation q± = qx ± q 2 to give 
normal modes 

q ± + (lo± V) 2 q ± = (45) 
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with eigenfrcquencics f2± = tu ± V, where we take the 
positive square root and assume \ V\ < u). As they should, 
these reproduce exactly the eigenenergies e± = e ± V of 
the quantum dimer problem. 

The momenta are obtained from the Hamilton equa- 
tions as 



which combine conveniently to give 



1 

P± = Pi ± P2 = q~9±- 



(46) 



(47) 



Thus in the eigenmodes we have the simple expression 
for the quantum amplitudes in terms of the classical am- 
plitudes 



so that in this case io s = uj — (K/2lo) = fi_. Thus, one 
takes two pendula of natural frequency u>. One then cou- 
ples them with strength K and simultaneously adjusts 
the lengths to give a new natural frequency u s . Then 
the oscillation amplitudes and velocities of this classical 
system reproduce exactly the complex time-dependent 
amplitudes of the quantum dimer wavefunction, with 
transition energy e = Huj and coupling matrix element 
V = K/(2uj). Clearly, as alternative one can leave the 
natural frequency unchanged as cu s and then, for given 
K infer the transition energy u> of the equivalent quan- 
tum system from uj s . Hence we have shown how to con- 
struct a pair of classical oscillators whose motions re- 
produce the entangled two-qubit quantum wavefunction 
timc-dcpcndencc. 



3. Dimer in RCA 



c± = z± 



Q± 



-Q± 



(48) 



Note that all the above results are invariant to a change 
of sign of V which only serves to flip the f2± eigenvalues. 

As example, in the beating mode the wavefunction is 
exactly reproduced by the classical amplitudes. The ini- 
tial conditions for beating are qi(0) = /?, 92(0) = 0, q\ = 
<72 = leading to amplitudes 



Qi/2(t) = (/?/2)[cos (Q+i) ± cos (0_i)] 



(49) 



Forming the velocities and from them the momenta, af- 
ter some algebra one finds P\/2{t) = — (/3/2)[sin (fi+i) ± 
sin (tl-i)] so that 



z 1/2 =(/3/2v / 2)[e 



-iQ+t 



±e 



(50) 



Choosing (3 = —y/2 and noting that e± = HQ±, these are 
exactly the quantum amplitudes c\ , C2 of Eq. (|41[) leading 
to the explicit beating forms of Eqs. (|42p . 

We now show how to simulate the quantum equation of 
motion using two coupled classical pendula. The classical 
equations of motion for two (mathematical, i.e linearized) 
pendula with oscillation angle 0, natural frequency ui and 
coupled by a spring with coupling strength K are 



£1 + W 2 (^i + Kfa = K(j) 2 
f>2 + U 2 4> 2 + K 4> 2 = K(f>i 



(51) 



Making again the identification V = K/{2uS) the equa- 
tions (|44[) which map exactly the Schrodingcr equation, 
can be written in the same form 



where 



q\ + uj s qi + Kqi = Kq 2 
q 2 + U) 2 s q 2 + Kq 2 = Kq 1 



K 2 /(4lu 2 



(52) 



(53) 



Previous works pointing out the equivalence of clas- 
sical oscillator motion and quantum time-development 
have not used the exact mapping but rather the simpler 
standard equations (|51|) or their equivalent and then in- 
voked the RCA. Next we investigate the accuracy of this 
approach in the simple dimer case. 

The eigenfrequencies of the exact mapping Eqs. (|52")) 



are readily calculated to be f2^_ 



and n 2 _ 



-2K. 



With K = 2coV this translates to the eigenfrequencies 
n± = oj ± V of Eqs. (HU and of course of the quantum 
problem. The eigenfrequencies of the standard Eqs. (|51[) 
are fi+ = u and i7_ = \/oj 2 + 2K involving the natural 
oscillator frequency. However in the realistic coupling 
approximation (RCA) we expand, as in section Til Bl 



= \JuP- + 2K kuj + (K/w) =uj + 2V 



(54) 



to give cigcnfrequcncy difference f2_ — 0+ = 2V, as in 
the exact mapping. This implies that, when the RCA 
is valid, the solutions to the standard Eqs. ([STj) will be 
a good approximation to the exact solution, up to an 
overall phase caused by the shift in absolute value of the 
eigenfrequencies. This explains the previously-observed 
close agreement of quantum solutions and standard (q- 
coupled) oscillator amplitudes. 

To assess the accuracy of the RCA we have calculated 
exact quantum (exactly equivalent to the solution of the 
p&q-coupled classical Eqs. (|4"4"|) ) and compared to solu- 
tions of the standard q-coupled classical Eqs. (f5"Tj) for the 
N = 2 dimer case with "beating" initial conditions. The 
results are shown in Figs. Q] and [2] In Fig. [1] the case 
w = 1, K = 0.01 is shown and one sees excellent agree- 
ment between exact, p&q- coupled (blue, solid) and stan- 
dard q-coupled (red, dashed) results over several periods 
of the transfer time (the time is given in units of K/(2u>) 
or equivalently V/H). Only a small shift in relative phase 
is perceptible. In the middle column the q-coupled results 
are multiplied by an overall phase factor cxp(i(f2_ — u))t) 
to compensate for the overall energy shift between clas- 
sical and quantum eigenvalues. The upper figures show 
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time in K/(2ui) time in K/{2uS) 



FIG. 1: Time dependence of the quantum or p&q-couplcd 
classical (blue) and q-coupled classical (red, dashed) motion. 
The left column shows details for short times and the right 
column is for longer times. The upper row shows the absolute 
value squared and the middle row the real part of the ampli- 
tudes. The bottom panel shows c*C2 and z*Z2. The time is 
given in units of K/(2ui) or equivalently V/h 



the absolute values squared of the classical q-coupled and 
quantum coefficients where this phase difference disap- 
pears. Similarly the lower panel shows the imaginary 
part of the off-diagonal density matrix elements c*C2 and 
z\z-i where the overall phase factor cancels also (the real 
part is zero). By contrast to the good agreement for 
K = 0.01, the case K = 0.1 is shown in Fig. [21 where 
one sees the beginning of the breakdown of the RCA in 
that, although still of the same shape, the phase differ- 
ence between quantum and classical q-coupled curves is 
becoming more pronounced. 

B. A circular array of N interacting monomers 

For N = 2 we have shown that by adjusting the nat- 
ural oscillator frequency to be wj one can readily con- 
struct a classical system to mimic exactly the quantum 
wavefunction. To illustrate the complexity of the equiv- 
alent classical system that begins to arise for quantum 
systems of larger Hilbert space dimension, we consider 
next the extension to N interacting quantum monomers. 
This is a model for Frenkel exciton energy transfer in 
molecular crystals [HI , on dye aggregates [1(1 [T3| or in 
the photosynthetic unit fl8j . In a previous paper [j| we 
have shown that in the RCA a set of q-coupled classical 
oscillators gives energy propagation characteristics indis- 
tinguishable from the exact quantum result. 



1. Quantum result 




±7T §7T §7T TV T 2?I " 

time in K/(2u) time in K/(2u) 

FIG. 2: Same as Fig.l except that the coupling strength is 
K = 0.1. 



For simplicity we take the monomers to be iden- 
tical with transition energy e. Furthermore we take 
only nearest-neighbour interaction into account and set 
V ra ,n±i = V. For circular boundary conditions the result 
of diagonalisation is standard [l6[ and gives 

E k = e + 2V cos k, (55) 

with k = (27r/7V)j, j = 0....N — 1. The transformation 
matrix elements are 

B kn = -^= exp(ikn). (56) 
v N 

If we consider the time development beginning with only 
monomer excited, i.e. c„(0) = 5 n o, then from the or- 
thogonality of normal modes one has 

Cn(t) = -jy ex P y— Jj i e ^J ^ cxp(^fcn) exp ^— -^(2V cos fc)i 

(57)' 

2. Exact mapping 

The exact mapping of this quantum system to the 
Hamilton equations gives the p&q-coupled Eqs. (|21l) in 
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the form 

q n + (uj 2 +2V 2 )q n = -2Vuj{q n+ i+q n -i)-V 2 (q n+2 +qn-2) 

(58) 

Note that, although the first-order Schrodinger equa- 
tions involve nearest-neighbour couplings, the equivalent 
second-order Newton equations explicitly contain next- 
nearest neighbour couplings. Substituting the trial solu- 
tion q n k = A kn exp (iflt + ikn) gives the eigenfrequency 
equation 

Q 2 - uj 2 - 2V 2 = AVuj cos k + 2V 2 cos 2k (59) 
with solutions 

n k = uj + 2V cos k (60) 

in agreement with the quantum result Eq. (|55[) . Hence 
for this iV-monomer case we see again how to construct 
a system of classical oscillators to reproduce exactly the 
quantum results. However, the couplings between oscil- 
lators is becoming more complicated. 

3. N-mer in RCA 

The simpler q-coupled classical equations arc those 
with nearest-neighbour coupling only, and were solved 
numerically in [5[ to compare with the quantum en- 
ergy transfer result. Here wc show analytically that in 
RCA one obtains the same eigenfrequencies and time- 
dependence as in quantum or p&q-coupled classical case. 
We consider a ring of identical coupled pendula. The 
coupling is between adjacent pendula only. The stan- 
dard equations of motion are, 

4> n + uj 2 <f> n + 2K(j) n = K((f> n+1 + n _i), (61) 

leading to the eigenmode equation 

[n 2 - (uj 2 + 2K)]B kn = -K(B k . n+1 + Sfc, n _i). (62) 

The eigenvalues are 

fl 2 k - (uj 2 + 2K) = -2K cos k, (63) 

which, in the RCA becomes 

fl k fn (uj + K/uj) - (K/uj) cos k. (64) 

Now we make the identification E k = Ml k ,e = hw and 
V/h= -K/(2uj) to give, 

E k = (e-2V) + 2Vcosfc. (65) 

Comparison to the quantum dispersion relation Eq. (|55j) 
shows an overall shift of eigenenergies by 2V. This is 
immaterial as it leads only to an overall phase factor. 
Hence, as was shown numerically in [§} for a particular 
set of initial conditions, the q-coupled classical oscillators 
describe to vey good approximation the exciton dynamics 
on an iV-monomcr chain. 



IV. CONCLUSIONS 

We have shown how the eigenenergies and eigenfunc- 
tions of an iV-level quantum system with real coupling 
matrix elements can be reproduced by a suitable array 
of coupled classical oscillators. Hence, by observing the 
classical motion, one can reconstruct the quantum me- 
chanical time-dependent wavefunction or density matrix. 
The difficulty of realisation of a classical system mimick- 
ing the quantum entanglement, is mainly in the condi- 
tion that a realistic system must reproduce all elements 
of the quantum couplings faithfully. This will be a task 
of increasing difficulty as the complexity of the quantum 
system increases. Nevertheless, we have shown explicitly 
that such a parallel is feasible for a circular or linear array 
of identical monomers. This is realised by a correspond- 
ing array of coupled pendula. Similarly, the excitation 
transfer on a molecular dimer, or similar coupled two- 
level quantum systems, can be faithfully simulated by the 
beating motion of a pair of classical oscillators, either me- 
chanical or electrical. Such experiments are carried out 
routinely in undergraduate physics laboratories. 

Realistically the multi-oscillator couplings necessary to 
construct p&q-coupled oscillators could be attained more 
easily with miniature LC oscillators than with mechan- 
ical ones. However, we have shown also that the com- 
plexity of p&q-coupled oscillators can be circumvented 
by using only q-coupled oscillators in the RCA. This ex- 
plains the success of previous comparisons of quantum 
time-dependent motion with that of physical q-coupled 
oscillator systems. 

Throughout we have neglected the effects of coupling 
of the quantum system or the classical oscillators to the 
environment. Of course in real systems such de-cohering 
or dissipating couplings are omnipresent. In the quan- 
tum case in the form of radiative or non-radiative decay 
of excited states or in eigenenergy shifts. In the classical 
case in the form of frictional forces for pendula or the 
equivalent resistance in the electrical case. At the phe- 
nomcnological level of simply assigning a complex energy 
to quantum levels or, equivalently, adding a term propor- 
tional to velocity to the classical equation of motion, it 
is easy to show that the quantum/classical equivalence 
is preserved (see Appendix [Cj . We have undertaken a 
more detailed study of this equivalence in the frame- 
work of standard theories of open quantum systems 0. 
Similarly, the correspondence between inducing radiative 
transitions in quantum systems and driving classical os- 
cillators with external fields is under investigation. 
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Appendix A: Commentary upon the RCA 

In the foregoing we have derived the RCA by approxi- 
mating the classical cigcnmodc frequencies. Further light 
can be shed on the nature of this approximation by look- 
ing at alternative derivations. A useful variation is to 
include explicitly the eigenvalue phases in the quantum 
time development of Eq. (|28[l i.e. to write 



(Al) 



*(t) = 5^a n (t)e- 



Forming the equivalent classical Hamiltonian as before 
we obtain 



U(t) = Y^H nm a* n {t)e 



ioj n t 



.We" 



(A2) 



where we have set tu n = e n /H. Note that if one proceeds 
to second quantization and elevates the coefficients a n to 
being annihilation operators, then this equation in just 
that of the quantum Hamiltonian in the Hciscnberg rep- 
resentation and the time-dependent phase factors are just 
the time propagators of the creation and annihilation op- 
erators. In particular, for fixed n and m the off-diagonal 
coupling terms are of the form, for H hermitian and real 



see that classical and quantum couplings will be the same 
in that the non-rotating wave terms can be neglected. 
A rather different view of the RCA is obtained from the 
Hamilton equations. The resulting q-coupled Newton 
equations are 



Qn + ^nQn — ~<^n H nm q m 



(A7) 



which arc identical in form to Eq. ([33]) for standard clas- 
sical oscillators. By contrast, for the p&q-coupling one 
has 



/ , H mra i q m ' , 



<jn+^n<ln — _ W n /J H nm q m — H nm 

m^n m^n rn' ^rn 

(A8) 

which differ from the Newton equations by the last term. 
However, again scaling as above we have 



K ^ 



(A9) 

Since in RCA we have ^ ;§> 1 we can neglect the second 
order term so that the p&q-couplcd hamiltonian gives 
the same second-order classical equations of motion as 
the q-couplcd hamiltonian. 



If one now sets 



I * 



—i(jjj n 



0*1 



(A3) 



a ne ^t = (q n (t)+ip n (t)) (A4) 

corresponding to classical amplitudes q n (t),p n {t) or 
equivalently Heisenberg operators, one obtains coupling 
elements 

T~Lnm = ~^H nm (q n q m + PnPm) (A5) 

which arc those of Eq. (|10|) . By contrast, if one restricts 
the coupling to terms q n q m as in the q-coupled case, one 
has 



(A6) 



y , — TJ [ * i(0J n -U} m )t , * -i{u) n -LJ m )t 

iL-nm — J ~ L nm Y n ^ c i u. m 1x77,0 



Then one sees that to obtain the hamiltonian of Eq. (|A3[) 
which maps exactly to the Schrodingcr equation, it is nec- 
essary to neglect the " non-rotating wave" terms involving 
the rapidly oscillating phase factors e ±, ( Un+u ' m ) t , which 
contribute weakly to transition probabilities. More ex- 
actly, let us take K to be the value of the largest of the 
elements H nm i.e. H nm — K"f nm where all 7 nm are less 
than unity. Then we can define a dimensionless time as 
t = Kt. The rotating wave factors are then of the order 
of exp(i(Aui / K)t) , where Aw is the mean frequency dif- 
ference and the non-rotating wave terms are of the order 
of exp(i(uj / K)t) , where u is the mean frequency. Now if 
bj/K 3> 1 and Atu/K w 1, which is the RCA, again we 



Appendix B: Coupled Electrical Oscillators 

As above, for the moment we ignore friction of resis- 
tance and consider an LC oscillator circuit. For a single 
oscillator the balance of e.m.f. from inductance and ca- 
pacitor is expressed by, 



4 + § = °- 



(Bl) 



where q is the charge and / is the current related by 



dt' 

Then we have the harmonic oscillator equation 



14+^ = 



or. 



w 2 q 







(B2) 



(B3) 



(B4) 



where the frequency is u> = \J\jLC. Note that, com- 
pared to a mechanical oscillator, L plays the role of mass, 
1/C that of coupling constant and —LI corresponds to 
momentum. 

Now we couple two identical LC oscillators by a ca- 
pacitor C K connected by leads in which current J = ^ 
flows. The equations of the coupled circuits are, 



d{h-J) qi 

dt C 

d(h + J) <12 

dt C 



(B5) 
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with the e.m.f. balance 

qi _ <12 _ Q_ 
C ~ C C K 

Using J = — one can eliminate J to obtain, 



(B6) 



lator equations (|33|) . ignoring the corresponding fluctua- 
tions, i.e. 



&n ~t~ 2X^ n x n -|- LO n X n ^ K nrn x ln . (^1) 



(1 + K)q\ +0J 2 qi - Kq 2 = 
(1 + K)q 2 +uj 2 q 2 - Kq\ = 0, 



(B7) 



where the dimensionless ratio K = Ck/C is dehned. 
These equations can be put also in the form of Eqs. (I5T1) . 



(1 + 2K)q\ + (1 + K)u 2 qi + Kuj 2 q 2 = 
(1 + 2K)q 2 + (1 + K)u 2 q 2 + Ku 2 qi = 0. 



(B8) 



Clearly when K — > the equations become those of un- 
coupled oscillators. Adding and subtracting either of 
these two sets of equations leads to uncoupled equations 
in the new variables q± = qi ± q 2l i.e. 



q+ 
q- 



n+q+ 

-.2 



=o 



n±q- =o 



(B9) 
(B10) 



which arc identical to Eqs. (|45|) except that now 0+ = 
u = y/T/LC and 0_ = u/y/{l + 2K). For coupled LC 
oscillators the RCA corresponds to K <C 1 to give, 



0_ 



u> — Kuj. 



(Bll) 



Comparison with Eq. (|54[) shows that in this case we 
make the identification 2V/H = Kuj. Then the coupled 
electrical oscillator equations giving q\ (t) , q 2 (t) are iden- 
tical to the q-coupled pendula equations and, as we have 
shown, in RCA reproduce the complex amplitudes of the 
quantum dimer. 



Then the coupled equations (|34|) become 



(nj -uj 2 n + 2in k T n )C kn = K nm C km . (C2) 



This equation is written in the form, 



ati. 1 ^ — ^ 
(f2fe— u) n -\-2i—- — ■ rr„)Cfe„ = — — ; y K nm Ckm- 

(ilk + UJn) (iifc + OJn) 

(C3) 

Again, in the RCA, we consider that the spread in eigen- 
frequencies VL^ and ui n is small compared to the mean 
natural frequency u> to approximate the above equations 
by 



(Oft U) n + iT n )Ckn — ~ — / K mn Ckm 

2uj t— 1 



(C4) 



If we assign each quantum level a width 7„ then the quan- 
tum coupled equations become, 



{E k - 6 n + i"fn)Bkn = ^ VnrnBk 



(C5) 



Appendix C: Inclusion of Damping 

To include damping phcnomenologically in the classi- 
cal case we add a velocity-dependent term to the oscil- 



which are identical in form to the classical equations. In 
Ref. Q it is shown, how pure dephasing can be realized 
in a system of coupled oscillators which reproduces the 
quantum results. 
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